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Abstract 

In this paper, we obtain optimal uniform lower tail estimates for the probability distribution of 
the properly scaled length of the longest up/right path of the last passage site percolation model 
considered by Johansson in [^. The estimates are used to prove a lower tail moderate deviation 
result for the model. The estimates also imply the convergence of moments, and also provide a 
verification of the universal scaling law relating the longitudinal and the transversal fluctuations of 
the model. 


1 Introduction 

In 1^, Johansson considered directed last passage site percolation on = {(m, n) : m,n G N} with 
geometric random variables. More precisely, for {i,j) G let w{i,j) be independent, identically 
distributed geometric random variables with 

nw{t,j) = k) = {i-e){ef, fc = o,i,2,..., (1.1) 

and 0 < t < 1. An up/right path tt from (1,1) to (M, N) is, by definition, a collection of sites {{ik,jk) ■ 
k = l,2,--- ,n},n := M + N-1 such that (b, ji) = (1,1), {in,jn) = {M,N) and (4+i, jfe+i)-(b, jfc) = 
(1,0) or (0,1). Let (1,1) /' (M, N) be the (finite) set of all such up/right paths from (1,1) to (M, N). 
Now define the maximal ‘length’, 

G(M, A) :=max{ ^ u;(f, j) : tt G (1,1) / (M, A^)}. (1.2) 

{i,j)G7r 
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Fix 7 > 1 and set M = [ 7 A^], the integer part of 7 -/V. The main result in is the following: for 
any a; G M, 


fGihN],N) - -^N 
lim P( J °° 

AT^OO \ boN^/^ 


< a; ) = F{a 


where 


1 


do — 


t((7 + l)t + 2v^)’ 


(1.3) 


(1.4) 


and 

(1-5) 

and where F(x) is the Tracy-Widom distribution j2^ for the largest eigenvalue of a random matrix 
chosen from the Gaussian unitary ensemble (GUE). In addition to ©, Johansson also proved large 
deviation results, 

hm ^ logP(G([ 7 fV], fV) < iV(- - y)) = -£(y), (1.6) 

N^oa JV^ Oq 

lim ^logP(G([ 7 iV],iV) > 7V(—-hy)) =-i(y), (1.7) 

N^oo i\ aQ 

for some explicit positive functions £(y) and i(y), y > 0. 

The result (|l.3D parallels an earlier result in [|| on the length of the longest increasing subsequence 
In (o') of a random permutation cr of fV letters. The main result in ||^ is the following: for any a: G K., 

where again F{x) is the Tracy-Widom distribution appearing in ( |1.3| ). The authors in also proved 
the convergence of moments: if XAr(cr) := and x is a random variable with distribution function 

F{x), then for m = 0,1, 2, • • •, 


hm E(x) 0 )=E(x" 
N—^oo 


(1.9) 


In earlier work other authors proved large deviation results for ijy, 

< v^(2 - y)) = -H{y), 


lim — 

N —»-oo N 
1 


N 


lin^-^logP(£Ar > '/N{2 + y)) = -I{y), 


( 1 . 10 ) 

( 1 . 11 ) 


for y > 0, where H{y),I{y) are certain explicit positive functions. The result (|1.10| ) is due to Deuschel 
and Zeitouni ||] and the result ( 1.11| ) is essentially due to Seppalainen |21 . 

In two recent papers, [Q [|^, the authors have considered in the moderate deviation regime. 
More precisely, for 0 < a < i, they showed that for y > 0 

iN<{2- yN-<^)Vi^ 


lim 

N^OO 


y^N 


l-3a 


1 

12 ’ 


( 1 . 12 ) 
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and 


lim 
N—*oo 


N > {2 + yN-^)VN) 


y3/2jv(l-3a)/2 

These moderate deviation results can be motivated by noting that 

d/12 


and 


F{x) 


1 - F{x) 


as x ^ — 00 , 




as X ^ + 00 . 


IhTTX^/^ 

Thus from (pT^), one anticipates that as iV ^ oo, 

<{2- yN-°‘)VN) = P(^Ar < 2Vn - {yN^/^-°‘)N^/^) 


l0gF(-2/iVl/3-«) ^ __y3^1-3a^ 


(1.13) 


(1.14) 


(1.15) 


(1.16) 


and 


P(^w > (2 + yN-°^)VN) = F{£n > 2'/N + {yN'^/^-°‘)N^^^) 

~ log(l - F{yN^^^-°‘)) ~ -ly3/2/v(i-3“)/2. 

3 


(1.17) 


Of course, when a = 0, we are in the large deviation regime, and when a = we are in the GUE central 
limit theorem regime. The above moderate deviation results can also be motivated by estimating the 
functions I{y) and F[{y) for the large-deviation regime. In ||^, the authors proved ( 1.13| ) by refining 
certain estimates in |Q| and using a careful summation argument. In [p^ , the authors utilized an 
analogous summation argument together with the estimate Lemma 6.3 (ii) in |Q|. 

Calculations similar to ( 1.16| ), ( 1.17 ), motivate the following moderate deviation results for G([ 77 V], TV): 
for 0 < a < |, 


lim 

N —»-oo 


and 


lim 

N —»-oo 


logP(G([ 7 iV],iV) < {^-yboN--)N) 

y3]\^2-3a 

logP(G([ 7 jV],TV) > {^+yboN--)N) 

y3/2 3oij2 


1 

12 ’ 


(1.18) 


(1.19) 


One of the principal goals in this paper is to prove (1.18). Relation (1.19) is slightly simpler and can 
also be approached using the techniques in this paper. We hope to return to this problem in the future. 


Relation (1.18) is a consequence of the following result. 

Theorem 1.1. Fix 0 < t < 1 and 70 > 1. Then there exist a (large) eonstant L > 0 and a (small) 
constant 5 > 0 , such that for large N, 

logP(G([ 7 TV],TV) < —TV - x6oTVi/3) = + oix'^N-^/^) + 0{\ogx) (1.20) 

Q-o 12 
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uniformly for all L < x < 5N'^/^ and 1 < 7 < 7 o- In particular, for the variables x, 7 in the same range, 

P(G([ 7 A^], A^) < —TV - xboN'^/^) < (1.21) 

Oo 

for some constant c > 0 . 

Setting X = in (1.20|), we immediately obtain (l.lSj) together with error estimates. 


Corollary 1.2 (Estimate for lower moderate deviation). For 0 < a < I and y > 0, as N ^ 00 , 


logP(G([ 7 iV],iV)<(i-y 6 oiV-“)iV) _ l ^ ^ , 

y3]\J2-3a 12 ”^^^ y3jg2-3a 


( 1 . 22 ) 


Theorem 1.1 can also be used for other applications. 


Corollary 1.3 (Convergence of moments). For 7 > 1, set Ojq := ^ o,nd let y be the 


random variable with distribution function F{x) as above. Then for m = 0,1,2, ■ 

hm E(0)O) = E(x™). 

AT— >CXD 


Proof. By Remark 2.5 of [^, ( 1.23| ) follows from the estimate (1.21). 


(1.23) 

□ 


In particular, setting to = 2, we see that the fluctuation Y^l/ar(G([ 7 TV],TV)) of G is of order 
where p — I. It is believed (see e.g., j^) that the transversal fluctuations of G have order TV^ where 
^ and 77 are related by a dimension-independent universal scaling law 2^ = 77 -|- 1. In other words, 
it is expected that in our case ■? = §• In (^, Johansson considered transversal fluctuations for the 
Poissonized version of the longest increasing subsequence problem and showed in that case that the 
scaling law 2^ = 77 -I- 1 is satisfied. By Q, 77 is again 1/3 and it follows therefore that the scaling law 
2^ = 7 + 1 is satisfied for this case. A key role in his analysis was again played by Lemma 6.3 (ii) of 1^. 
This Poissonized problem can be viewed as a continuum version of the above site percolation problem 
and in Remark 1.2 of Johansson notes that the scaling law 2^ = 77 -|- 1 for the site percolation 


problem would follow from an estimate of type (1.21) above. The modifications in the argument in |13| 
that are needed for the site percolation problem are detailed in jlj]. We thus have 

Corollary 1.4 (Transversal fluctuations). For any 7 > 1, the above coefficients 77 and ^ for longi¬ 
tudinal and transversal fluctuations of the site percolation model obey the scaling law 


2 ^ = 77 + !. 


(1.24) 


In order to prove Corollary 1^, IT, weaker bounds than ( 1.21 ) suffice. Indeed, using an observation 
of Harold Widom (see in particular Lemma 2), it is possible to prove the bound 


P(G([7A^], A^) < —N - xboN^/^) < 

tto 


(1.25) 


for x, 7 in the range of Theorem 1.1, for some constant d > 0. As opposed to the proof of (1.21), which 
requires a steepest-descent Riemann-Hilbert analysis (see below), the proof of (1.25|) uses only classical 
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steepest-descent methods. A key role in is played by a beautiful conjecture of Widom for the 
spectral properties of a class of singular integral operators (see @). This conjecture can be verified in 
our case, as in the case considered considered by Widom in Q , by using an elegant formula of Borodin 
and Okounkov (see identity (4.9)). The method in is itself motivated by earlier calculations in 
[||. The estimate ( 1.25| ) is enough to prove Corollary O, but does not suffice to prove Corollary 

O. 


Remark 1. The bound in Lemma 6.3 (ii) of [Q was also used by Seppalainen |2^ to control fluctuations 
for the “stick process” introduced in [Q. In |^, Seppalainen also mentioned that a similar result 
could be obtained for a certain continuous-time totally asymmetric simple exclusion process, provided 
the appropriate analogue of Lemma 6.3 (ii) could be established. The same should be true for a discrete¬ 
time version of this process. The above estimate ( 1.25 ), and of course also the stronger estimate ( 1.2l| ), 
then suffices to control the fluctuations as in the stick process. 

Remark 2. Our results are given for G{M, N) where M = [yA^] > N, but as the statistics of G{M, N) 
are clearly the same as for G{N,M), it is immediate that our results, suitably scaled, also apply to 
G([ 7 A^],fV) for 0 < 7 < 1. 


As indicated above, the proof of Theorem 1.1 is based on the steepest-descent method for Riemann- 
Hilbert problems (RHP’s) introduced by Deift and Zhou §, and further developed in [Q. The method 
has been used to solve a wide variety of asymptotic problems in pure and applied mathematics (see, for 
example, |^, ^ and the references therein). The steepest-descent calculations in this paper are closely 
related to the calculations in IJ and particularly [Q|. Our analysis is based on the algebraic formula 


(2.7) below, which relates P(G(M, A) < n) to the solution Y to an associated RHP on the unit circle 


S = {|z| = 1} (see (^)). It follows then that our problem reduces to the asymptotic analysis of a 
RHP with large oscillatory parameters. The steepest-descent method in Q was introduced precisely for 
this purpose. The key step in the method is to identify the leading order asymptotics for the solution 
of the RHP and this is done following i I and by introducing a so-called g-function with 


certain specific properties on an appropriate contour Pi U r 2 (see Proposition 4.1 below). Using g, one 
transforms the RHP for Y as follows: U := where ^ is a specific constant to be 

determined and is the Pauli matrix (T 3 = (J ^ simple calculation shows that U solves the RHP 


U{z) is analytic in z G C \ E, and continous up to the boundary, 
U+iz) = U-iz)\ ,, , I, zeE, 


(1.26) 


0 


=fe(s+-s-) 


where W is given ( ^.17 ), and g± denote the boundary values of g. In addition, one requires g{z) = 
logz -I- 0(z“^) as z ^ 00 , so that the RHP for U is normalized at infinity, 

U{z) = I + 0{1/z), as z —> 00 . (1-27) 

The choice of the properties of g mentioned above is made precisely such that the leading contribution 
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to the RHP (1.26) is immediate. Further information on the steepest-descent method can be found, for 
example, in @], ||, ||, §. 

In [Q, |j, the RHP’s are given on the real line R and the analogues of Pi, r 2 are subintervals 

of K: in [Q, the RHP is given on the unit circle S = {\z\ = 1} and the analogues of Pi, r 2 are again 
subintervals of S. The main new technical feature of the RHP in this paper is that Pi and r 2 cannot 
be chosen as subintervals of the original contour E, and the central problem is to discover the shape 
of Fi, F 2 . The situation is similar to the problem confronted by Kamvissis, McLaughlin and Miller in 
[ p^ , where the authors considered the semi-classical limit of the solution of the Cauchy problem for 
the focusing nonlinear Schrodinger equation. Motivated by the calculations in jl^, we show that the 
construction of Fi, F 2 is equivalent to the problem of determining the global structure of the trajectories 
Q{z){dz)‘^ > 0 and orthogonal trajectories Q{z){dz)‘^ < 0 of a particular quadratic differential Q{z){dz)‘^ 
(see (3.53) below). 

The outline of the paper is as follows. In Section H, we derive the basic algebraic formula (2.7) 
relating P(G([ 77 V], TV) < n) and the RHP (^), and state our basic asymptotic estimate. Proposition 
2.2 , for Y2i{0;k). In Section ||, which is the heart of the paper, we construct Ei and S 2 using the 
theory of quadratic differentials and verify the desired properties of h = g'. In Section |^, the constant 
£ mentioned above is defined (see (4.2)) and the desired properties of g = f^h are verified (Proposition 
4.1). In Section |^, we use the 5 -function to analyze the RHP and eventually give the proof of Proposition 


2.2. Finally, in Section g, we use the estimate in Proposition 2.2 together with a careful summation 


argument as in ll^ to prove the main result Theorem 1.1 
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2 Basic relations and formulae 

For M,N> 1, let 

( 2 . 1 ) 

Set 

^(z) = (l + tz)^(l + l)^, ( 2 . 2 ) 

z 
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and consider the n x n Toeplitz determinant 


— Dji :— det((/^j_/j;)o<j,fc<n5 


(2.3) 


where (pj is the Fourier coefficient of (p: 


ipj := 



z ^p>{^) 


dz 

2Triz 


(2.4) 


Here and below the integration contour \z\ = 1 is assumed to be oriented in the counter-clockwise 
direction. Let G{M, N) be the maximal length introduced in the Introduction. From earlier result of 
Gessel an Johansson ||^, Baik and Rains Q extracted the relation 

P(G(M, N)<n) = -i—(2.5) 

^M,N 


which plays the basic role in our analysis. 

Let S be the unit circle \z\ = 1 in the complex plane, oriented counter-clockwise and let Y(z) = 
y(z; k) = (Fy (z; fc))i<ij <2 be the solution to the following 2x2 matrix RHP: 


F(z) is analytic in z G C \ E, and continous up to the boundary, 
1 z-'^^iyz) 


Y+{z) = Y_{z) 


z G E, 


.0 


1 


^y(z)^-fc<T3 = 7 + o(l/z), asz^oo. 


( 2 . 6 ) 


where = (o-?i) is the standard third Pauli matrix, and Y+{z), (resp., F_(z)), z G E, are the 
boundary values of Y{z') as z' ^ z from the inside (resp., outside) of the circle. 


Lemma 2.1. The solution Y to the above RHP (|2.6|) exists and is unique. Moreover, 


P(G(M, iV) < n) = n (-F2 i(0 ; k + 1)). (2.7) 

k—n 

Proof. We will construct the solution Y explicitly using computations similar to ||l^ . First note that 
from the equality (|2.5|), Dn{p>) 7 ^ 0 for n > 0 since the probability P(G(M, iV) = 0) = V{w{i,j) = 0,1 < 
Vi < M, 1 < Vj < N) = (1 — , and hence P(G(M, TV) < n) > P(G(M, iV) = 0) > 0 for n > 0. 

Consider for k > 0 the polynomials of degree k 


/ To 

T-1 ■ 

T-k \ 


( To 

Ti 

Tk \ 

Ti 

To 

T-k+1 

> K{z)-=^Aet 
T>k 

T-i 

To 

Tk-1 

Tk-l 

Tk-2 ■ 

T-1 


T-k+l 

T-k+2 

Ti 

\ 1 

z 

z'^ / 


V z'^ 

zfe -1 . 

■ 1 / 
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A direct check shows that TTfc and satisfy the following orthogonality conditions: 


l\z\=l 

f dz 

/ z~^Trl{z)(p{z)-^— = NkSjo 0<j<k. 

J\z\ = l 


z ^Trk{z)ip{z] 


dz 

2mz 


= NkSjk 0 < j < fc, 


where 


Nk = 


Dk+i 

Dk 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


Let {Ch){z) = /s ^ C\S, denote the Cauchy transform of h. Let {C±h){z) = limzi^z{Ch){z') 

where z' approaches z from the ± side respectively, denote its boundary values as in the Introduction. 
We claim that 


Y{z;k) = 


T^k{z) 


C'((-) Wk){z) 


-N^\K_,{z) -iV,-_\C((.)-V7rLi)(z) 


( 2 . 12 ) 


is a solution to (2.6). The analyticity of y in C \ S is clear, while the jump condition follows from the 


relation (7+ — C- = 1. The asymptotic condition follows from the orthogonality (2.9), (pTQ). On the 


other hand, the uniqueness of the solution to the RHP (2.6) is standard (cf. for example. Lemma 4.1 
of 1^). Hence ( 2.12 ) is the unique solution to the RHP (2.6). 


For the proof of (2.7), first note that by taking n oo in (2.5), 

lim Dn{^) = Zm,n- (2-13) 

n—^■oo 

(This can also be seen directly from the Szego strong limit theorem for Toeplitz determinants.) Thus 


we have, using (2.11), 


OO oo 

P(G(M, IV) < n) = n ^ = n 

k—n k—n 


Finally, from (2.8) and ( 2.12| ), we observe that T 2 i( 0 ; k) = which completes the proof. 


Remark 3. From (2.7), 


~^ 2 i( 0 ; fc + 1 ) — 


and hence we have — 121 ( 0 ; k) > 0 for fc > 1. 


P(G(M,7V) < k) 
P(G(M,iV) < fc + 1)’ 


(2.14) 


□ 


(2.15) 


There are three parameters in the RHP (p(^): M, N, k. We regard t, 0 < t < 1, as a fixed number 
throughout this paper. For convenience we introduce the following notation: 


7 := 


M 

iV’ 


N 


a := 


(2.16) 





















Instead of M, N, k, we now regard 7 , a, k as our parameters in the RHP ( |2.6D . Using this notation, the 
jump condition for Y takes the form 
/1 g-fcw^(z)\ 

y+(z) = F_ (z) ( ^ j, IU(z) :=- 7 alog(l + tz) - alog(l + t/z) + logz, (2-17) 

where logic is defined to be analytic in C \ (— 00 , 0], and logic = log |ic| for ic > 0. We are interested 
in the asymptotics of U as k 00 , while 7 and a remain in appropriate bounded regions. The main 
technical result we are going to prove in the rest of this paper is the following: 


Proposition 2.2. Set 


1 

t((7 +l)t + 2 ^)' 

Fix 0 < t < 1 and 70 > 1- There are Lq, Jq, fco > 0 such that for 7 , a, k satisfying 

flo + < a < (1 + (5o)aoi 

k > ko and 1 < 7 < 70 , we have 

log(-y' 2 i( 0 ; A:)) = -C 2 k{a - ao)^ + 0{k\a - ao|^) + 0{\a - ao|) + 0{— -r), 

k\a — ao| 


where 

^ t^(t + t7 +2^7)3^ 

4(l + tV7)2(t + ^)2- 

In particular, for the variables in the same range, 

log(-l2i(0; k)) < -cok{a - aof 


(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


for some constant Cq > 0. 


3 The /i-function 


As noted in the Introduction, we seek a ^-function and associated contours Pi Ur 2 . Rather than 
analyzing g directly, we first seek an h-function, h' = g. We start from an ansatz, to be verified a 
posteriori, that Pi, r 2 has the shape indicated in Figure with the endpoints f, Recall from (2.17) 


W{z) = — 7 alog(l + tz) — alog(l + t/z) + logz, 
and recall from (2.18|) 


W'{z) = - 


7 a 


z + t 


-1 


z + t 


+ 


n T 1 


(3.1) 


t(( 7 +l)t + 2 ^)' 

Notation from the theory of quadratic differentials (see e.g., Given a meromorphic function /(z) 

and a simple oriented curve C lying outside the zeros and poles of /, the notation f{z){dz)^ > 0 
means that /(z(t))(^)^ is real and positive for all t G {a,b) where z{t), t G (a, 6 ) is the arc length 
parameterization of C. Similarly, f{z)(dz)'^ < 0 means that /(z(t))(^)^ is real and negative. 
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Proposition 3.1. Fix 0<t<l, 7 > 1 , a> cq. Then there exist (cf. figure [J below) 

• a point ^ with Im(^) > 0 , 

• a simple open curve Fi connecting ^ and oriented from ^ to ^ such that (i) it does not intersect 


(— 00 , 0 ] and (ii) Fi is symmetric with reflection about the real line, 


• a simple open curve F 2 connecting ^ and oriented from f to ^ such that (i) it does not intersect 


(— 00 ,—1/t] U [—t, 00 ) and (ii) F 2 is symmetric with respect to reflection about the real line, 
• a function h{z) analytic m C \ Fi and continuous up to the boundary 
such that the following properties are satisfied: 

(a) h+{z) + h-{z) = W'{z) for z G Fi. 

(h) h{z) = j + 0(z~^)) as z ^ 00 . 

(c) i{h+{z) — h-{z))dz > 0 for z G Fi. 

(d) {2h{z) — W'{z))dz < 0 for z G F 2 ft C+ and (2/i(z) — W'{z))dz > 0 for z G F 2 n C_. 

In addition, we have the following properties: 

(i) The function h has the form 



(3.3) 


where R{z) 


{z — £)){z — 0) is defined to be analytic in C \ Fi and R{z) ~ z as z ^ 00 . 


(ii) (Endpoint condition) 



(3.4) 


(Hi) Set 


$(z) := (1 + ja) 



(3.5) 


where 


-i?( 0 )(l + 7 a)' 


(3.6) 


Then 


2/i(z) - IF'(z) = 4>(z), zGC\Fi. 


(3.7) 
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In particular, 


h+{z) - h-{z) = <I>+{z) 
2h{z) - W'{z) = $(0), 


z G Fi, 

(3.8) 

z G r2. 

(3.9) 


(iv) The curve r2 intersects the real axis at —zq. 



Figure 1: Sketch of the contours Fi and F2. 


Idea of proof : 

Suppose that the curve Fi is known. Let 


R{z) = J{z - Oiz - 0, 


(3.10) 


which is defined to be analytic in C \ Fi, and R{z) ~ 2: as z ^ 00. From the Plemelj formula, it is easy 
to check that 


h{z) = 


R{z) 


W'is) 


-ds, 


2Tri Jn R+{s){s - z) 
satisfies condition (a). By a residue calculation, h can be written as 

7a ( 


h{z) = \w'{z) + \R{z) 


u -j- 1 


(z + t ^)R{—t {z-\-t)R{—t) zi?(0) 

For this h to satisfy (b), the following two conditions are necessary and sufficient, 


W'{s) 


ds = 0, 


1 


_ sW'js) 

2-Ki R+(s) 


Vri R+is) 

or equivalently by residue calculations, 

7a a a + 1 

R{-t) ^ i?(0) ’ R{-t) 

As we will see, these conditions determine the endpoint 


ds = — 1, 


7a 


at 


= -(1 + 70). 


(3.11) 


(3.12) 


(3.13) 


(3.14) 
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Now for z e Ti, from |U^) , 

i{h+{z) — h-{z)) = iR^{z) ( -^ + -r—-—r-■ (3.15) 


Substituting R{—t and R{—t) in terms of i?(0) using the endpoint conditions ( 3.14 ) we obtain, after 
some algebra, for 2 G Fi, 


i{h+{z) — h-{z)) = i{l + 7a)i?+(z) 


{z + Zo) 


z{z + t){z + t“^) ’ 


where 


Hence we have 


zo = 


Qj \ 


Q{z) := {i{h+ - hyy = -(1 + 7a) 


-i?(0)(l + 7a)' 

2{z - £.){z-0{z + zy 


(3.16) 


(3.17) 


(3.18) 


Z2(z + t)2(z + t-l)2 ’ 

which is meromorphic in C. We then use the theory of the quadratic differentials to find the trajectories 
for Q{z){dzy > 0 and this leads us to the determination of the contour Fi for which condition (c) is 
satisfied. 

The contour F2 for which condition (d) is satisfied turns out to be obtained by finding the so-called 
orthogonal trajectories corresponding to Q{z){dzy < 0. Clearly if a trajectory and an orthogonal 
trajectory meet at a point z in the plane where Q{z) is analytic and nonzero, they do so at right angles 
to each other. 

The rest of this section consists of a proof of the above Proposition. 

3.1 The endpoint ^ 

In this subsection, we are going to prove that there is a unique Im{^) > 0, for which the following 


two conditions (cf. ( 3.13 )) are satisfied, 

Jt R+{s) 


ds = 0, 


1 f sW'{s) 
2TTi Jy R+{s) 


ds = —1, 


(3.19) 


where F is any simple oriented curve connecting ^ and oriented from ^ to ^ such that (i) it does not 
intersect (—oo,0] (ii) it is symmetric under reflection about the real line, and 


R{z) = \J {z - - i), 


(3.20) 


which is defined to be analytic in C \ F with R{z) ^ z as z —> 00. 

Remark A priori we should look for a pair of unrelated points ^i, ^2 such that ( 3.19| ) is satisfied for 
any contour F connecting them as above. It turns out, however, that it is sufficient to look for and 
^2 in the form Ci = and ^2 = The reason for this symmetry lies in the form of the equations (3.19|). 
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Indeed, both W'{s) and sbb'(s) are real analytic functions, and in each integral the path of integration 
can be doubled along the “minus” side of the branch cut for R{s) and then deformed into a closed loop 
containing and ^2 and the branch cut F connecting them but no singularities of W'{s)] this loop is 
otherwise arbitrary. If we take the loop to be symmetric with respect to reflection in the real axis, then 
it is easy to see that the only way for both integrals to be purely imaginary as required by (3.19) is for 
i?(s) itself to be a real analytic function, which forces ^2 = ■Ci- 


As in (3.14), these conditions become 

7 a a a + 1 


7 a 


at 


i?(-t-i) R{-t) i?(0)’ ti?(-t-i) R{-t) 


= -(1 + 70 ). 


Set 


l?l =-^(0): x:=\t +^\=-R{-t ), y:=\t + ^\ =-R{-t). 


Then conditions (3.21) have the form 


7 a 

X 


a 

y 


7 a at 

- - 1 -= 1 

tx y 


■ 7 a. 


Set 


ri := 


t{a + 1 ) 

1 + 7 a ' 


r2 := 


u 1 
t(l + 7 a) 


The conditions (3.23) are now, after simple algebra, equivalent to 


X 

a(l — t^) 

y 


= t(l+ 7 a)(l - — 


= t(l+ 7a)(—- 1). 


Also from the definitions, r, x, y satisfy the relation 


2 ^ 

r = 1 + 


l-f2 


Inserting x,y of ( tl.25| ), ( 3.26| ) into ( 3.27 ), we obtain an equation for r: 


r 


1 , (i-t> 

2 


1 


-1 = 0 . 


Since x,y > 0, from (3.25), 


(1 + 70 )^ \t^(r 2 — r)^ {r — ri)'^ ^ 

, we must have ri < r < r 2 . Thus we seek x, y, r satisfying 

X > 0 , y > 0 , n < r < r 2 . 


(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 


Lemma 3.2. For each fixed 0<t<l, 7>1, a> oq, there is a unique solution r to (3.28) satisfying 


ri < r < r 2 - 
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Proof. Set 


H{r) := ^ + 


1 , 


2'i„2 


1 


- 1 . 


(3.30) 


(l+ 7 a)^ \t^(r 2 —r )2 {r — ri^ ^ 

Clearly, H{r) —oo as r J, ri, and ff(r) ^ +(X) as r t 'f 2 - Thus there is ri < < r 2 satisfying 

H{rc) = 0. We want to show that such an is unique. By direct calculation, for ri < r < r 2 , 




1 


tY 


H{r)+^H'{r) (i + ^o )3 - r)3 ' (r - ri)3 

The minimum of this function on (ri, r 2 ) is obtained at 

_ ri + r2y/yt 


- 1 . 


1 + V7t ’ 


and for ri < r < r 2 . 


H{r) + -H'ir) > H{r.) + = 


a(l + Vlt) 

(l + a){l-t 2 ) 


2 \ 2 


- 1 . 


But since a > Og, 


H{r) + ^H'{r) > 


ao(l + yYt) 


2 \ 2 


-1 = 0 . 


(3.31) 


(3.32) 


(3.33) 


(3.34) 


,(1 + ao)(l - 

Therefore, if H(rc) = 0 for ri < Tc < ^ 2 , we must have H'frc) > 0. A simple calculus argument then 
proves the uniqueness of the solution. □ 


Thus if we define y by ( [1.25| ), ( 3.26 ), we have obtained the unique solution (r, a;, y) to the equations 
( 3.25 ), ( |3.26 ), ( 3.27 ) subject to ( 3.29 ). Now we need to prove that the {r,x,y) defined in this way 
determines > 0, uniquely from (3.22). In order for Im{£) > 0, to satisfy (3.22), we must 

have ^ = re*® for some 0 < 6 * < tt satisfying 


COS0 = 


X — — t 


2 y2 _ _ ^2 


2rt 


-1 


2rt 


(3.35) 


Conversely, if there exists 9 G ( 0 , tt) satisfying ( |3.35D , then ^ := re“ and ^ = re “ are the desired 
endpoints. However, the second inequality follows from (3.27) and so it is sufficient to prove that for 


(3.36) 


{x,y,r) satisfying ( |3.25D , ( p.26| ), (|3.27D , ( p.29[ ), we have the relation 


x^ — r'^ — t ^ 

2rt- 


In order to prove (3.36), we first prove the following Lemma. 


Lemma 3.3. For each 0<t<l, 7 > 1 , a> og, the solution {r,x,y) to ( 3.25 ), ( 3.26 ), ( 3.27 ), subject 

(3.37) 


to (3.29) satisfies 


X + y > t ^ —t. 
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Proof. From ( 3.25| ), ( |3.26 ), we have 

x + y 


t ^ — t 1 + 70 


7 — 1 + 


ri7 


r2 


r — ri r 2 — r 


(3.38) 


The minimum of the right-hand side, regarded as a function in r, is again obtained at r = r*, where r* 


is defined in (3.32), and hence, by evaluating the minimum, we obtain 


x + y ^ a{t + V7)^ 


But since a > oq, we have 


-t (1-1- 7a)(l - t^) ’ 


x + y ^ apit + ^f 


t ^ -t (1 -I- 7ao)(l - f^) 


(3.39) 

(3.40) 

□ 


Now we prove (3.36). 


Lemma 3.4. For each 0 < t < 1 , 7 > 1 , a > oq, the solution {r,x,y) to (|3.25D , (|3.26D , (|3.27| ), subject 

(3.41) 


to (3.2£) satisfies 


-1 < 


x^ — r^ — t ^ 


2rt 


-1 


< 1 . 


Proof. Suppose that 


x^ — r"^ — t ^ 


2rt 


-1 


> 1 . 


Then x"^ > {r + t ^)^, and thus from ( |3.27 ), 

^ (^2 _ _ ^ 2 ) > ^ 2 (^ ^ ^- 1)2 ^ (^2 _ _ ^ 2 ) ^ ^ ^^2 ^ 

Hence we have 


X > r + t 


-1 


y>r + t. 


(3.42) 


(3.43) 


(3.44) 


Inserting (3.44) into (3.25), (3.26), we have 






> t{l + ja){^ - 1 ) 


(3.45) 


, ,_i - > i(l + 7a)(l-—), , , 

r-|-t r ' r +1 

We multiply the first inequality by r(r -|-1“^), and multiply the second inequality by r(r -|- 1). Then by 
adding the two inequalities, we obtain, after some algebra, 0 > 2 , which is a contradiction. 

Now suppose that 


2 2 4.-2 

X — r —t 

2rt~^ 


< - 1 . 


This implies that (te)^ < (1 — rt)'^. Since r < r 2 < \ for 7 > 1, we have 

tx < \ — rt. 


(3.46) 


(3.47) 
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and from ( |3.27 ), 


y 2 = ^ 2^2 (-^2 _ _ ^ 2 ) < (1 _ + (r^ - 1)(1 - _ i) 2 ^ (-3 

Hence we have y < |r — t|. We distinguish two cases r >t and r < t. For the first case when r > t, 
from (3.47) we have 

X + y < — r) + (r — t) = — t, 

which contradicts Lemma |3.3| . For the second case when r < t, 

y <t — r. 

Inserting ( 3.47 ), ( 3.50| ) into ( 3.25 ), (|3.26D , we obtain 


(3.49) 


(3.50) 


07(1 — t^) t(a + 1 ) a(l — t^) 0 + 1 , 

<1+70--^- -7 - -< -t(l+7a). 


(3.51) 


1 — rt ' r ’ t — r r 

Multiply the first inequality by (1 —rt)r, multiply the second inequality by (t—r)r, and add the resulting 
two inequalities. Then after some algebra, we hnd o < 0, which is a contradiction. This proves the 
lemma. □ 


It now follows from the preceding discussion that there is a unique solution /m(^) > 0, to ( 3.21 ), 
or equivalently ( 3.19| ). 


3.2 The contour Fi 


As in (|3.6| ), set 


zo ■■= 


0+1 


(3.52) 


-i?( 0 )(l + 7 o)’ 

where i?(0) is given in ( 3.22| ). We emphasize that zq is uniquely determined by the endpoint and 
is independent of the curve F in subsection |3.1| , as long as F does not intersect (—oo,0], Note that 
t < Zo < from the condition n < r := —i?(0) < r 2 of ( 3.29|) . 

Now define 


Q{z) := -(1 + 70 )^ 


{z - Qjz - 0{z + Zq)^ 

z^iz+tnz+t-^r 


(3.53) 


and we consider the quadratic differential Q{z){dz)'^. In this subsection, we are interested in the 
trajectories Q{z){dz)‘^ > 0. In the next subsection we consider the orthogonal trajectories Q{z){dz)‘^ < 
0. A general reference for quadratic differentials is Chapter 8 of [^. Note, in particular, that if two 
trajectories (or two orthogonal trajectories) meet at a point z where < 5 ( 2 :) is analytic and nonzero, then 
they must be identical. The quadratic differential Q{dz)‘^ has double poles at and 00 , a 

double zero at —zq and simple zeros at We first consider the local structure of the trajectories near 
the poles and the zeros. 
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3.2.1 Local structure 

• Near 0: Q{dz)'^ ~ ^^{dzY, where c = (1 + 7 a)|^|zo > 0. Thus we want trajectories z = z{t) 

satisfying ^ ci for some real constant ci. The solution is given by z ^ and hence 

the trajectories near z = 0 are circular. 

• Near —t: Q{dz)^ ~ (dz)'^, where c = (^I > 0. As in the z = 0 case, the 

trajectories are circular. 

• Near Qidz)"^ ~ (dz)'^, where c = —±£oi ^ q_ Again, the trajectories 

are circular. 

• Near 00 , Q{dzY ^ -^{dz^ = -^{dw)^, where c = 1 + 7 a > 0, ic = A Thus once again the 
trajectories are circular. 

• Near —zq: Q{dz)‘^ ~ —c^(z + zo)^(dz)^, where c = > 0. Thus we seek trajectories 

z = z{t) satisfying ic{z + zq)^ ~ Ci for some real constant ci. The solution with z(0) = —Zq is 
given by (z + Zq)^ ^ —i(2ci/c)t, thus we have arg{z + Zq) = j + fc = 0,1, 2, 3. Hence there 
are 4 trajectories starting from —zq, all making an angle with the real line. 

• Near Q{dz)‘^ ~ c(z — ^)(dz)^ for some c G C. Thus we seek trajectories z = z{t) satisfying 

Cy^z — ^ Cl for some real constant ci. The solution with z(0) = ^ is given by (z — ^ 

(ci/c)t, and hence arg{z — C) = C 2 + fc = 0,1, 2. Thus there are 3 trajectories starting from 

making an angle ^ between themselves. 

• Near Q{dz)‘^ ~ c(z — ^){dzY for some c G C. Again, there are 3 trajectories starting from 
making an angle ^ between themselves. 

The local structure of the trajectories of Q{dz)^ near the poles and the zeros are summarized in 
Figure ||. 


3.2.2 Global structure 


Note that as Q{z) = Q(z), if {z{t) : a < t < (3} is a trajectory of Q{dz)^ > 0, then {z(t) : a < t < P} 
is also a trajectory. Also from ( ^.18 ), we see that Q{z) < 0 for z G K \ {—zo,0, —t, hence all 

trajectories that cross the real axis do so at 7r/2, and if {z{t) : a < t < P} is a trajectory that satisfies 
Im{z{t)) > 0 for a < t < /3, and Im{z{P)) = 0, Re{z{P)) —zo, —t, 0, then {z(2/3 — t) : P < t < 

2/3 — a} gives a smooth continuation of {z{t)} into Im{z) < 0. 

We need the following lemma. 


Lemma 3.5. Let^,^ be as in subsection 3.1. LetT be a simple curve with endpoints which does not 


intersect (— oo,0] and is symmetric under reflection about the real axis. Choose the branch of \JQ(z) 

i(l+7a) 


to be analytic in C \ F, and satisfy Q(z) 


as z 00 . For any real number x satisfying 
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Figure 2: Local structure of the trajectories of Q{dz)^ near the poles and the zeros. 

< X < —t, let C be a smooth curve in C-|- with endpoints ^ and x, oriented from ^ to x, which 
does not intersect F (see Figure^. Then we have 

Re [ y/Qdz = 0. (3.54) 

Jc 





Figure 3: The contours C and C. 

Proof. Let C be the closure of CUC, oriented from ( to Hence C" is a curve which has endpoints 
intersects the real axis at —t~^ < x < —t, and satisfies C" = C' (see Figure ||). Let C* = {z : z G C}, 
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oriented from ^ to x. Using ^/Q{^ = —'jQiz) and the realness of x, we have 

y/Q{z)dz + [ y^Q{z)dz 

Jc* 


Re [ \jQ{z)dz = ^ / \/Q{z)dz + [ y^Q{z)dz 

Jc ^ IJC Jc 


\/ Q{z)dz + / -JOi^dz 


uc 


l-c* 


Uc 

1 


(3.55) 


^Q{z)dz. 


Hence we want to prove that the last integral is 0. 
Set (cf. dlTil) , (PI) ) 




W'{s) 


-ds 


r R+{s){s - z) 

- 2 ^ (^) + V(2 + t-i)i?(-t-i) + (z + t)i?(-<) “ Ip)) j 

i?(z)(z + Zo) 


(3.56) 




where the second equality follows from a residue calculation, while the third equality follows from 
the endpoint conditions (3.21). Thus from the definition ( 3.53| ) of Q and the choice of VQj we have 
iQ = i{2h — W). Now 


[ 'JQ{z)dz = [ {2h — W')dz = [ dz [ — rds— [ W'{s)ds 

Jc Jc Jc ttz J^R+{s){s-z) Jc 


1 r W'{s) 
ni Jr i?+(s) 

By a residue calculation, for s € T, 


ds [ —^-^dz — [ W\s)ds. 
Jc s-z Jc 


(3.57) 


5 <£>,, = ta 1 

r—^■oo 2 


R{z) 


Ic S- z 


s — z 


>|=r 


dz — 7rii?_|_(s) = TTl 




— s — R+{s) 


(3.58) 


Hence using the endpoints conditions ( p.l9 ), we have 

1 - . - «+(.))* - W'{s)ds 

= 2Tri — J W'{s)ds — J W'{s)ds 
= 2-Ki — J W'{s)ds, 


(3.59) 


where T', the closure of T U C, encloses —t, 0, but not —t and is oriented counter-clockwise. But a 
direct calculation shows Jp, W'{s)ds = 27ri, and we obtain the lemma. 


□ 


Remark 5. The authors are indebted to Nick Ercolani who suggested that a formula such as ( 3.54 ) 
should be true. 

The following general results for the trajectories of quadratic differentials are given in Lemmas 8.3, 
8.4 of 1^, and are used in several places in the arguments that follow. 
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Figure 4: Local structure of the trajectories of Q{dz)^ near —zq. 


Lemma 3.6. Let Q(z){dz)‘^ be a quadratic differential in a simply connected domain G. 

(i). If there is at most one pole of Q{z) in G and this pole is simple, then there is no closed Jordan 
curve in G consisting only of trajectories (or orthogonal trajectories) and their endpoints. 

(ii). Suppose that Q{z) has no poles in G and let T be a trajectory (or an orthogonal trajectory). Then 
in both directions, F ends at a zero of Q or converges to dG. 


Remark 6. At various points in the argument that follows, we will use Lemma 3.6 (i) in a slightly 
extended form. For example, we will want to consider Jordan curves consisting of trajectories of the 
form in the first picture of Figure ^ for which the hypotheses of the Lemma are not fully satisfied. 
However, if we, for example, make a change of variables z ^ = z^~^, 0 < e < 1, then the figure 

takes the form as the second picture in Figure ^ and as the change of variable takes trajectories to 
trajectories, it is easy to verify that the hypotheses of Lemma |3.6| (i) are now satisfied. We will use this 
extended form of the Lemma without further comment below. 




Figure 5: Unfolding of a “closed” trajectory meeting itself at a point. 


Denote the four trajectories emerging from —zo by 1,2, 3,4 as shown in Figure ^ First, consider 
the trajectory emerging from —zq along the ray 1. Since there are no poles in C+, which is a simply 
connected region, this trajectory must either go to the zero of Q in C+, or escape from C+. Suppose 
that the ray 1 does not escape from C+. Then from the definition of a trajectory, '/Qdz G ]R\{0}, and 


hence the ray 1 (and similarly the ray 2) can not go to ^ due to Lemma 3.5. So the ray 1 must exit from 
C+. Now it can not exit through —zq because then by the local structure of the trajectories, the ray 1 
comes back to zq through ray 2, and the interior of the loop has no poles of Q, contradicting Lemma 
3.6 (i). Also, again by the local structure of the trajectories, the ray 1 can not exit through —t~^,—t, 0 
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or oo. Now there are five possibilities: the ray 1 exits from C+ through (—oo, —t ^), (—t zq), 

(—zo,—i), 0), or (0,oo). We examine each case. 

(i) The ray 1 can not exit through (—zq), for if it does, the trajectory can be continued by 
complex conjugation as remarked before, and the interior of the loop contains no poles of Q{dz)'^, 
contradicting Lemma |3.6| (i). 

(ii) By the same argument, the ray can not exit through {—zq, —t). 


(hi) Suppose the ray 1 exits at zi G (—t, 0). Then extending the ray by conjugation, we obtain the 
closed loop (see the first Picture in Figure Now the trajectory along ray 2 from —zq also can 
not go to ^ and hence must exit C+. Also the ray 2 can not cross the trajectory emerging from 



Figure 6 : Ray 1, case (hi). 


the ray 1. Hence it must cross the real axis between (—zo, —t) or (—t, zi). The first case can not 
hold as the closed loop obtained by extending the ray 2 by conjugation has no poles inside. In 
the latter case, the simply connected region formed by the two loops has no poles, which is again 
a contradiction (see the second picture in Figure ^). Thus ray 1 can not exit through (—1,0). 

(iv) Suppose that the ray 1 exits at Z 2 S (0, oo). Then by arguing as in case (hi), the trajectory along 
the ray 2 from —zq must exit C+ as some point Z 3 € (—t, 0). We denote the open regions in C+ 



Figure 7: Ray 1, case (iv). 

divided by the trajectories emerging from —zq by I, II, III as shown in Figure]^. Now the zero ^ 
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of Q lies in one of the regions I, II, III. (Note that ^ can not be on the trajectory emerging from 
—zo due to Lemma 3 . 5 .) 


(iv-1) If C G III, then as before, each of the three trajectories emerging from ^ exit through 
{—zo,—t) or {—t,zz). Hence at least two of the trajectories exit through the same interval. 
Then the simply connected region bounded by these two trajectories has no pole, which is 
again a contradiction. 

(iv-2) If ^ € II, we reach a contradiction by a similar argument. 

(iv- 3 ) If ^ G J, the three trajectories emerging from ^ must exit, one through (—oo,—one 
through Z4 G {—t~^,—zo), and one through (z2,oo). But this in turn gives a contradiction, 
because Re ^/Qdz = 0 by Lemma 

Hence in each case, we have a contradiction. 

Thus we conclude that the trajectory along the ray 1 emerging from —zq exits C+ through Z5 G 

Now we consider the trajectory along the ray 2 emerging from —zq. As before, it must exit C+ 
through either of (—00,25), 0) or (0, 00). 

(i) If it exits through (—00,25), the two trajectories from the rays 1 and 2 , extended as in the 
remark above to C+, form a simply connected region which does not contain poles, this is again 
a contradiction. 

(ii) Suppose that the trajectory along the ray 2 exits through {—zo,—t). Then the (C_-extended) 
loop of the trajectory 2 contains no poles, and again we have a contradiction. 



Figure 8: Ray 2 , case (hi). 

(hi) Suppose that the trajectory along the ray 2 exits through zg G ( 0 , 00). Then ^ lies either in the 
region I, H or HI as of Figure ||. If ^ G I, at least two of the trajectories from ^ exit together 
through (—00,25) or (20,00). This yields a contradiction as in (iv-I) above. The case ^ G III 
leads to a similar contradiction. If ^ G II, then the three trajectories exit, one through (—20, —t), 
one through {—t, 0 ), and one through ( 0 , zg). But this gives in turn a contradiction by Lemma 
as in the case (iv- 3 ) above. 

Therefore the trajectory along the ray 2 emerging from — zg must exit C_|_ through 27 G {—t, 0 ). 
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Figure 9: Rays 1 and 2. 


Thus ^ lies either in regions I, II or III of Figure But by a now familiar argument as above, ^ 
can not lie in II or III. Hence ^ S /. Then the three trajectories emerging from ^ exit at some points 
zs,zg,zio where zg € (—oo,zg),zg € (zy,0),zio G (0,oo). This shows that the global structure of the 
trajectories of Q{dz)'^ is given in Figure [Toj 



Figure 10: Global structure of the trajectories of Q{dz)^. 


The above considerations prove, in particular, that there is a trajectory emerging from ^ and ending 
at 5 which crosses the real axis to the right of 0. We take Fi to be this curve, oriented from ^ to 
Define h by ( 3.56| ) with the choice of the contour F = Fi: let R{z) = \J{z — ^){z — C) be analytic in 
C \ Fi, R{z) ~ z as z ^ oo. Thus 


h{z) = 


R{z) r W'{s) 

2TTi in R+{s){s - z) 


ds 


= Iw'iz) + ii?(z) 


70 


R(z)(z + zo) 


(z + t)R{-t) 


t + I \ 


(3.60) 




Condition (a) of Proposition |t[ h+{z) + h_{z) = IF'(z) for z e Fi, now follows by properties of the 
Cauchy operator, and the condition (b) by the endpoint condition ( tl.igj ). 

Now we consider condition (c). Define Q{z) = z(I + 70 ) analytic in 
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C \ r 1 , and ^/Q 


z(l+7a) 


as 2 ^ oo. Then 


i{h+{z) - h-{z)) = \/Qiz)^, zeTi. 


(3.61) 


Along the trajectory Ti, we must have y/Qdz € R \ {0}. Let zio G (0,oo) be the point at which Ti 
crosses the real axis. As R{z) = {z — ^){z — > 0 for z G (zio,oo), it follows that for all upward- 
oriented trajectories that cross the real axis at points x G (zio,oo), ^/Qdz\z=x < 0. It particular, for 
Ti, VQ-dz\zi^ < 0, and hence \jQ{z)_dz < 0 for all z G Ti. Now since VQ+ = —VQ- on Ti, we have 

^/Q^dzGR+, onTi, (3.62) 


and hence from (3.61), 


z(/i+(z) — h_(z))dz > 0, z G Ti, 


(3.63) 


which proves condition (c). 


3.3 The contour r 2 

Again we choose the branch of \/Q so that ^/Q is analytic in C \ Ti and -y/Q ~ di+ 70 ) as z ^ 00 . 
Now we consider the orthogonal trajectories Q{z){dz)‘^ < 0. As before, the local structure is easy to 
determine. We summarize the local structure of the orthogonal trajectories near the finite poles and the 
zeros of Q{dz)^ in Figure ^ Near 00 , any straight rays emerging from 00 are orthogonal trajectories. 
We note that the real axis is an orthogonal trajectory. Hence the orthogonal trajectories can cross the 





Figure 11: Local structure of the orthogonal trajectories near the finite poles and zeros of Q{dz)^. 
real axis only at 0 , —t, Zq, —t~^ or 00. 

Now we consider the global structure. Again by the symmetry Q{z) = Q(z), the orthogonal tra¬ 
jectories are symmetric under reflection about the real axis. There are three orthogonal trajectories. 
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denoted by 1, 2, 3, emerging from each of which bisects the angle between two trajectories emerging 
from ^ (see Figure |I^. Now we show that the orthogonal trajectory 1 can not cross either of the two 



Figure 12: Orthogonal trajectories and trajectories emerging from 

adjacent trajectories emerging from If so, there is a closed loop of the form shown the Figure 
consisting of a part of a trajectory and a part of an orthogonal trajectory, and by analyticity of '/Q, 
the integral of ^/Qdz is zero around the loop. But the integral of ^/Qdz along the trajectory is in 

2 


-^-•- 

Z9 0 

Figure 13: Orthogonal trajectories emerging from 

]R\ {0}, and along the orthogonal trajectory, the integral of ^/Qdz along the orthogonal trajectory is in 
iR\ {0}. Hence the sum can not be zero, which is a contradiction. Therefore the orthogonal trajectory 
1 must exit C+ between zg and zig. But then from the local structure, it must exit at 0. Similarly, the 
orthogonal trajectory 2 must go to cx). The orthogonal trajectory 3, by a similar argument, must exit 
C+ either through —or —t. Suppose it exits through —t. By the local structure, the orthogonal 
trajectory approaches along an angle; indeed it is easy to show that z{s) ~ as s ^ 00 for 

some 0 < (/) < TT, c > 0. Also by the definition of an orthogonal trajectory, we have 

s) _ 

\/ Q{z)dz = 0 (3.64) 

along the orthogonal trajectory 3. For s large, but fixed, write z(s) := —t + ee*® for e > 0, 0 < 0 < tt, 
and consider the curve Cq from ^ to z{s) along the orthogonal trajectory 3. Let Cg be the curve 
{—t + ee*^ ■ d < P < tt}, oriented from z{s) to —t — e. Thus Cq U is a curve from ^ to a point in 


lim Re 
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{—t Then by Lemma we have Re^/Qdz = 0. Hence we have 


lim Re / \/Qdz = lim Re / \/Qdz — Re / \/Qdz 
JCo JCqUCs JCs 

= — lim Re / \/Qdz 

JCs 

rJ + * 

’ Jc. 


(3.65) 


= — lim 


z{z + t~^) z + t 




This contradicts (3.64), and hence the trajectory 3 cannot exit C+ through —t. Similar argument shows 
that it cannot exit C+ through —t~^. Hence the trajectory 3 exits C+ through —zq. Thus we have 
obtained the global structure for the orthogonal trajectories of Q{dz)^ emerging from the zeros 
—zq as shown by the solid curves in the Figure M- The dotted curves in the Figure 0 denote the 
trajectories already displayed in Figure M. 



Figure 14: Global structure for the orthogonal trajectories and the trajectories emerging from the zeros 
I and -zo- 


In particular, we have shown that there is an orthogonal trajectory (more precisely, a union of two 
orthogonal trajectories, one in C+ and the other, its conjugate, in C_, meeting at the point —zq which is 
a zero of Q{dz)‘^) emerging from ^ and ending at ^ and which crosses the real axis at —zq. We denote the 
curve by r 2 , and take the orientation from ^ to Since r 2 is a union of two orthogonal trajectories, we 
have y/Qdz S i]R\{0} for z G r2\{0} and an explicit computation using + 

shows that 


y/Qdz G zM-, 
^/Qdz G zK+, 


on F 2 n C+ 


(3.66) 

(3.67) 
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on r2 n c_ 









(3.68) 

(3.69) 


with the orientation from ^ to From ( |3.60| ), i{2h — W') = ^/Q. Thus we have 

i2h-W')dz<0, zGr 2 nC+, 
i2h-W')dz>0, zGFsnC-. 


This proves condition (d). The reader will observe that the remaining conditions and formulae in 
Proposition 3.1 have been proved en route in this section, and this completes the proof of the Proposition. 

These facts can be illustrated by numerical computations of trajectories and orthogonal trajectories 
associated with the quadratic differential Q{z) (dz)^. The (orthogonal) trajectories can be obtained 
simply with a Runge-Kutta scheme. For (orthogonal) trajectories that emerge from zeros or poles of 
Q(z), which amount to singularities of the vector field in the complex plane, the only additional difficulty 
is to determine the initial direction, which is not unique. But the possible initial directions are easily 
determined by the sort of local analysis that has already been presented above. An example of the 
results of such a calculation is presented in Figure p^. 


All Trajectories: a = 4.00000, 7 = 2.00000, t=0.707107 



Figure 15: The trajectories and orthogonal trajectories of the quadratic differential Q{z) {dz^ emerging 
from the points z = ^, z = and z = —zq- The parameter values are a = 4, 7 = 2, and t = l/-\/2. The 
trajectories where Q{z) {dz^ > 0 are shown with thick curves and the orthogonal trajectories where 
Q{z) {dzY < 0 are shown with thin curves. 

A computer program for generating numerical approximations to the contours Fi and r 2 is of course 
useful because it allows one to explore/illustrate the dependence of the contours on the parameters a, 
7 , and t. As an example. Figure ^ illustrates the deformation of the contours as a is varied while 7 
and t are held fixed. Here, we see clearly what happens as a is decreased to qq, its minimum value 
which depends on 7 and t. Namely, the two endpoints converge to a common point on the negative real 
axis, and at the same time, the contour Fi closes without collapsing, while r 2 disappears. Similarly, 
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Con;ourb T. ard r.: a=0,150000, 7=2,00000, t=0.707107 Contours r. ord r.: o=0,500000, 7=2.00000, t=0.707107 Contours r, ord r,: o=1.00000, 7=2,00000, t=0.707107 



Figure 16: The dependence of the contours Fi and r 2 on the parameter a. As a tends to ag from above, 
the endpoints ^ and ^ coalesce on the negative real axis. 

as a increases without bound, the opposite situation prevails, with the endpoints ^ and ^ coalescing on 
the positive real axis, while r 2 closes without collapsing and Fi disappears. It is also possible to see in 
these pictures that the closed curve Fi U r2 deforms somewhat throughout this process; the endpoints 
^ and ^ do not simply slide along a fixed closed curve as a is varied. 


4 The ^r-function 


Proposition 4.1. Fix 0 <t<l,7>l,a>ao. Let Ti, r 2 be as in the Proposition 3.1. Let Inz 
denote the principal branch of logarithm, Inz G M for z > 0, and set 

9{z) ■= [ {h{s) --)ds + lnz, z G C \ (Fi U (-oo,pi]), (4.1) 

d oo ^ 

where pi > 0 is the intersection point o/Fi and R, and the integral is taken over a curve from oo to z 
which does not intersect Fi U {—oo,pi\. Let 4>(z) be as in (3.5) of Proposition 3.1. Also set 

£-.= 2g{i)-W{i). (4.2) 


Then g and £ satisfy the following properties: 

(1) g{z) is well-defined and analytic in <C\ (Fi U (—oo,pi]), and is analytic inC\ Fi. 

(2) = z(l + 0(l/z)) as z ^ oo. 

(3) e3+(A+9-{A-w{z)-e = i for z e Fi. 
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(4) 9-{z) = g'i'i+( 2 ) joj. 2 . G Fi, where = exp{/^^ $(s)(is}, z G C\riUR_ with the 

integral taken along any curve not intersecting Fi U R-. 

(5) = g'i' 2 (z) Jqj. 2 g F 2 , where '^ 2 {z) = $(s)ds with the integral taken along F 2 . 

Remark 7. It follows from ( |3.3| ) that h{z) = h{z) for z G C \ Fi, and hence from ( |4.l| ), we see that 
g{z) = g{z) for z G C \ (Fi U {- 00 ,pi)). 

Proof. Since h is analytic in C \ Fi and continuous up to the boundary, we have g given by @) is 
well-defined and analytic in C \ Fi U (—oo,pi]), and continuous up to the boundary. Now let C be a 
simple closed curve, oriented counter-clockwise enclosing Fi. Flsing the formula for h, 


h{z)dz = [ ^^dz [ 


bF'(s) 


-ds = 


j_ r f ^ 

ic Jc Jri -R+(s)(s - z)““ 27ri7ri^+(s) Jcs-z 

Using the residue at infinity, we then have 


dz. 


1 


27ri !r s — z 


R{z). f + f 

-dz = s — 


2 ’ 


and 


/c=/r. 


(4.3) 


(4.4) 


(4.5) 


from the endpoint conditions (^). Now from the definition of g, for z G (—oo,pi), g+(z) — g-(z) = 
Jg, h{s)ds = —2'Ki, where g± (z) = limej^o g{z ± ei). Therefore is analytic in C \ Fi and continuous 
up to the boundary. This proves property (1). 

Since h{z) = - -\-0{z~‘^) as z ^ 00 , we have g{z) = lnz-|-0(z“^) as z ^ 00 , which proves property 

( 2 ). 

For z G Fi n C+, 

9±i.z)=9{0+[ {h±{s) --)ds+ \nz-Inf = g{f) + f h±{s)ds, (4.6) 

s 

where the integral from ^ to z is taken along Fi. Hence from h+{z) + h-{z) = lU'(z), z G Fi, we have, 
using (4.2), 

g+iz)+g.iz) = 2gif) + Wis)ds = W{z) + £, (4.7) 

for z G Fi n C+. For z G Fi n C_, 

9±{z)=9(i) + J_ h±{s)ds = g(f) + J_ h±{s)ds j h±{s)ds, (4.8) 

where the integrals are again taken along Fi. But from ( |4.5| ) (recall that h is analytic in C \ Fi), 
Jj{h^{s) — h_{s))ds = 27ri and g{f) — g{f) = jJh_{s)ds. Hence for z G Fi n C_, 


g±{z) = 9 {f) — J h-{s)ds J h±{s)ds J h±{s)ds 


(4.9) 
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and hence 


g+{z) + g-{z) = 2g{^) + / {h+{s) — h-{s))ds + / W'{s)ds = W{z) + £ + 2TTi. 


(4.10) 


Therefore, since e®+, e®- and W are continuous for z S Ti, 
z € Ti, which verifies property (3). 


From (3.8), (4.£), (4.9), and the above relation /|(h+(s) 


we have = 1 for all 

- h-(s))ds = 2Tri, we have for z G Fi, 


„g+(z)-g_(z) ^ expjy (h+(s) - /l_(s))ds| = + 


(4.11) 


and hence the property (4) follows if we prove that exp{ ^{s)ds}, z G C \ Fi U IR_ does not depend 
on the choice of the integration path. For this purpose, it is enough to prove that 


$(z)dz G 27riZ 


(4.12) 


for any simple closed contour C which encloses Fi and does not intersect (—oo,0]. From (0) and the 
fact that W'{z) is analytic away from —1,0, we have 


Using 


/ $(z)(iz = / (2/i(z) — W'{z))dz = 2, [ h{z)dz. 

Jc Jc Jc 

for h, the above integral is equal to 


(4.13) 


i?(z) 


dz 


lU'(s) 


rds = — 
ni 


»-■'(*),, f m,,f *)*. ( 4 . 14 ) 


Jc Jti R+is){s - z) ni R+{s) JqS-z R+i-s) 

But this is equal to 47ri from the endpoint condition (3.4), and so ( 4.12| ) is established. 
For z G r 2 n C+, as in (|4.6|) 


9{z) = giO + J h{s)ds 


(4.15) 


and 


2g(z) - VF(z) -£ = 2 h{s)ds - W{z) + 1U(C) = ^ {2h{s) - W'{s))ds. (4.16) 

For z G r2 n C_, as in (^), 

g{z) = g{£,) + _ h{s)ds = g{^) — _ /i_(s)ds+ _ h+{s) + / h{s)ds = g{^) + / h{s)ds + 2ni, 

(4.17) 

and hence for z G r2, 

.. ^ ^ (4.18) 


^2giz)-w{.)-e ^ g^pj J\2h{s) - VF'(s))ds|, 


and property (5) follows from (3.9). This completes the proof of Proposition 4.1 


□ 
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5 RHP analysis 


Set 


r = riur2, (5.1) 

oriented counter-clockwise. It is a simple closed curve which has 0 and —t inside and —t~^ outside. 


Since the jump matrix ) for Y in (2.6) is analytic in C\ {0}, we can deform the contour 

S for Y to r, as follows. Set 


Y{z) for z inside both F and S, and for z outside both F and E, 
Y{z) = F(z)t4'(z) for z inside F and outside S, 

Y{z)Vy^{z) for z outside F and inside S. 


(5.2) 


Then Y is analytic in C \ F and continuous up to the boundaries, satisfies I + (z) = Y-Vy(z) for z S F, 
and Y{z)z~^'^^ = I + 0{z~^) as z ^ oo. 

Now (see the Introduction) we define 


M(z) = 

Then from Proposition |4.l| (1), (2), M satisfies 
M(z) is analytic in z S C \ F, 
M+(z) = M_(z) 


g-k{g+{z)-g-iz)} ^k{g+{z)+g-{z)-W{z)-e) 


0 gHg+i^)-a-i^)) 

M{z) = I + 0{l/z), as z ^ oo. 

From the Proposition |4.l| (3), (4), the jump matrix V for M is now 

/g-feifi+lz) Y \ 


z G F, 


V{z) = 


0 


ofc'I'i+(z) 


z G Fi 


and from the Proposition (1), (5), we have 

/1 pfc’I'2(z) 

t/(z) = 


z G F 2 . 


^0 I ) 

For the jump matrix on z G Fi, note that ^^+( 2 :) = —4'(z)_ and 

/g-fc>i'i+(z) 1 \ / I o\ / 0 1 


V{z) = 


0 g-fe'i'i-(d 


1 O' 

g-fe'&i_(z) i] \_1 qI L-fe'&l+(z) J 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


Clearly 'I'i+ has an analytic continuation to the (-l-)-side of the contour Fi. Now for z G Fi, it is 
easy to see that i?e('I'i^(z)) = 0, and hence from the Proposition 3.1 (c) and (3.8), the derivative 


31 









^4'i+(z) = $+( 2 ;) along the contour Fi satisfies /to(^'I'i+(z)) < 0, z G Fi. Thus the Cauchy- 
Riemann condition for the analyticity implies that Re{'^i{z)) > 0 for 2 : on the (+)-side of Fi and close 
to the contour. Therefore we can take a contour F^^^ with endpoints ^ for which i?e('I>(z)) > 0 
for 2 ; G int(r^^^). Similarly, 'I'i_ has an analytic continuation to the (—)-side of Fi and its real part 
is positive for z on the (+)-side of Fi close to the contour Fi and we take a contour F^ ^ for which 
i?e(d>(z)) > 0 on its interior. We take the orientation of F^, j = 1,2 to be from ^ to See Figure E 
for the general shape of the contours F^^^ F^^\ 



Figure 17: The contours F^^^ and F^^^ 


Define M{z) to be M{z) for z in the region bounded by F 2 and F^^^ and also in the unbounded 
region. For the region bounded by Fi and fJ^\ define 


M = M 


1 


O' 


J’ 


( 2 ) 

and for the region bounded by Fi and F^ , define 


M = M 


1 


O' 


g-fe'&FD I ■ 


Set F' = Fi U F^^^ U Fp^ U F 2 . Then M satisfies the new RHP 

M(z) is analytic in z G C \ F', 
M+{z) = M_{z)V{z) zGT', 
M(z) = / + 0(l/z), as z ^ 00 , 


(5.8) 


(5.9) 


(5.10) 
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where the jump matrix V is 


V(z) = 



z G Fi, 


z G urf\ 


Z G Ft. 


(5.11) 


Now we take the limit fc —*■ oo with 7 > 1 in a compact set and ag < a < a*, for some a*. From 
the signature of i?e(d>i(z)) on Fp\ j = 1, 2, we see that the jump matrix V ^ I. For z G F 2 , we have 
'i' 2 {z) < 0 from the Proposition lA and Indeed ^> 2 ( 2 ) is decreasing as z follows from ^ to —zq 
along F 2 , and then increasing as z follows from —zq to ^ along F 2 . On the other hand, d>2(0 = 0 
'^' 2(0 = = 0 , as $(s)ds = ^{s)ds = —$(s)cfs, and the negativity of 'I' 2 (z) on F 2 follows. 


Hence as k 00 , V 


I on F 2 . Therefore we have V 
0 1 

V°° = l\ -1 0 


V°° where 


z G Fi, 


(5.12) 


zGF^^^UFf^UF2. 


Let M°° be the solution to the RHP with the jump matrix V°° and normalized at infinity. The solution 
is given by 

/ §+£: 


M°°iz) = 


0-0- 

2i 


where 


/3(z) = 




0 - 0 - 

2i 

0 + 0 - 


1/4 


i" 2 l( 0 )' 


,-fe(s( 0 )-£) 


M“(0), 


(5.13) 


(5.14) 


which is defined to be analytic C \ Fi and satisfies /3(z) ^ 1 as z ^ 00 . We expect that M ^ M°° as 
k 00 , and hence by tracking the algebraic transformations Y ^ M ^ M, we expect that 


(5.15) 


for 7 > 1 in a compact set and ag < a < a*. 
Indeed we have: 


Proposition 5.1. Let 1 < 7 < 71 for any fixed 71 > 1. There are Li, i5i > 0 such that for 


have 


for sufficiently large k. 


UQ + ^2/3 — ® — (1 + <^l)®0) 


T2i(0)e-'=(s(°)-^) = + 0( 


k\a — Qg 


■)), 


(5.16) 


(5.17) 
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The convergence V —>■ V°° is not uniform on T' and this considerably complicates the analysis. As 
in ^ in order to obtain the above error bound, we need to introduce local parametrices for the solution 
of the RHP near each of the endpoints ^ and As in |]p, a suitable local parametrix near each endpoint 
can be obtained in terms of Airy functions. Also since a is not fixed, but is allowed to approach uq, we 
need to vary the magnitude of the parametrix according to the size of a — oq. A similar situation arises 
in Lemma 6.2 (ii) of |^. The proof of the above Proposition is parallel to that of Lemma 6.2 (ii), [Q ( 7 , 
in Q play the same role as a, k in this paper, respectively), and we do not repeat the argument here. 


Lemma 5.2. We have A := ^(0 + iO) — £ € K. In particular, ^ ^ > 0. 

Proof. For a; G M \ {pi}, by Remark |^, g{x + iO) = g{x — iO), and hence by Proposition 4.1 (1), 
gS(x+*o) ^ gg(x-*o) ^ and so is real. In particular, it follows that 

Im{g{x + zO)) G Ztt, and hence by continuity, Im{g{x + zO)) is constant for x < pi. From (iD and 
the proof of Proposition O (3), £ = g+{C) + 9-{f,) — bF(C) = 9+{Pi + *0) + 9-iPi) ~ W{pi). Hence 
gs(o)-f _ gg(o+iO)-g+(pi+iO)g-g-(pi)+w(pi)_ above, lm{g{0 + zO)) = Im{g+{pi + zO)). Hence 

g{0 + zO) — g+{pi + Oz) G K. Clearly g-{pi) and W{pi) are also real, and this proves the lemma. □ 

If we set ^ 0 < 9c < tt, then we can check /3(0) = and Mff{0) = —13{0)~^) = 

— sin ^. Hence the above proposition yields that 


->^ 2 i( 0 ) = e 


fcA 


siny(l + 0 (-r-i 

2 ^ k\a — ao\ 


))■ 


(5.18) 


Note that — 121 ( 0 ) is indeed real and positive from Lemma |2.l| . This is consistent with Lemma 5.2. 
Now we compute Let 


^ ^ lf“l n 

a = —— = | 4 | cos 02 . 


(5.19) 


From the formula (O) for h, one can check directly that an anti-derivative of 2{h{z) — 1) is 
- 7 alog(z + t~^)- olog(z -I- £) -I- (a - 1) log(z) 


7 a 

X 


G “h 1 


R{z) - ( - — (a + t ^) - -(a -I- £) -f 
X y 


n -f 1 


— 7 a log 


z -I- R{z) +t ^ — X 
z + R{z) + + X 


— a log 


z - 1 - i?(z) + t — y 
z + R{z) + t + y 


-I- (a -I- 1 ) log 


a ) log(z - a -I- i?(z)) 


z -k i?(z) - r 
z-ki?(z)-kr 


(5.20) 


where the logarithms are taken to be the analytic in C \ (—oo,0] and logz = log|z| for z > 0, and 
r = —i?(0), X = —R{—t~^), y = —R{—f) as in ( ^.22| ). It is straightforward, but tedious, to check that 
(5.20) is analytic in C\ (Fi U (—oo,pi]) as in the definition of g in (^t|). Using the endpoint conditions 
(3.4), or( 3.23| ), ( 5.20 ) is equal to 


- 7 olog(z -k £ ^) - olog(z -k £) -k (a - 1 ) log(z) -k (1 -k 7 a) log(z - a -k Riz)) 


— 7 a log 


z -k R{z) +t ^ : 

Z “k R{z) t ^ “k ^ 


— a log 


z -k i?(z) + t-y 
z -k i?(z) + t + y 


-k (a -k 1 ) log 


z -k R{z) — r 
z + R{z) + r 


(5.21) 
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Evaluating the asymptotics as z ^ oo, we see that 

2 g{z) = - 7 alog(z + t~^) - alog(z + t) + (o + 1 ) log(z) 
+ (1 + 7 a) log((z - a + R{z))/2) - 7 alog^ 

' z + R{z) + t-y 


— a log 


z + R{z) + t + y 


+ (a + 1 ) log 


z + i?(z) + — X 

z + R{z) + + X 

z + i?(z) - r 


R{z) 


and 


2Re{g{0 + lO)) =(7 - l)alogt + (1 + 7 a) log((Q; + r)/2) 

—r + t — y 


— 7 a log 


—r + t ^ — X 


—r + + X 

Also from (p. 22 ) and (4.2), we have 


— a log 


—r-\-t + y 


+ (a + 1 ) log 


2r 


1 + a/r 


t = 7 alogt + (1 + 7 a) log(^ - a )/2 

+ f^ + t-y 


— 7 a log 


^ + t-^ + X 


- alog 


^ + t + y 


+ (a + 1 ) log 


^ + r 


By (B) , 


cos 9c = 


X — r — t 


2 rt 


-1 


y2 _ _ ^2 

2 rt 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


Thus using (5.23) and (5.24), we can express A in terms of r,x,y. After some algebra, we find 
A = — 70 log t + (2 + a + 7 a) log 2 + (4 + a) logr — i log(r + — x) 


1 


1 


(5.26) 


- 0(1 + 27 a) log(r + t + a;) - - log(r + t-y)- -{1 + 2 a) log(r + t + y). 


To emphasize the dependence on o, we write A = A(a), etc. 

Lemma 5.3. Fix 4 < 72 < 00 . Then there exists (52 > 0 such that for oo < a < (1 + 62)00 and 
1 < 7 < 72 , we have 


A(a) = -C 2 (a - ao)^(l + 0{\a - ao|) 


where 


_ t^(t + t 7 +277)^77 
4(l + tV7)2(i + V7)2’ 

and the order term 0 {\a — ao|) is uniform for a and 7 as above. 


(5.27) 


(5.28) 


For the proof, we need the following lemma, which considers the case when a = oq. This case is 


specifically excluded from Lemma 3.2 
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Lemma 5.4. For fixed 0 < t < 1, 7 > 1, and a = uq = > there is a unique solution r to 

( 3.28 ) satisfying ri <r <r 2 , given by 


r = r(ao) = rg := 


l + iVl 


i + Vl 

The function a 1-^ r(a) is smooth for all a > ag and 

3(7- I)(t7 + t + 2^)2t2 


r'{ao) = — 


4(1+ ty7)(t+ 77)3 


(5.29) 


(5.30) 


Proof. When a = ag, we have ri = trg and r2 = 


t'o- 


Since t < rg < j, rg satisfies ri < rg < r2. 


It is then a direct calculation to check that rg satisfies the equation ( 3.28 ). Now we want to show the 
uniqueness of the solution. Let H be as in ( ^.30| ) in the proof of Lemma |3.2| . Now when a = ag, the 
value r* of (3.32) is rg. Thus we have 


H(r} + ^iJ'(r) > R(rg) + ^H'(rg) = 0, 


(5.31) 


for ri < r < r 2 , and the inequality is strict for r ^ rg. Hence if there is a zero Vc ^ rg, it should satisfy 
H'(rc) ^ 0. On the other hcincl, by direct Ccilciilntion, we hcive 

H'(ro) = = 0, (r,) = ^ 

and hence H is also increasing at the zero r = rg. As in Lemma |3.2| , there is no other zero rc 7 ^0 in 

{ri,r2). 

Now consider H = H{a, r). By direct calculations, we find 


H{ag,rg) = Hr{ag,rg) = Ha{ag,rg) = Hrr{ag,rg) = 0 


(5.33) 


and 

4 (^ + y/l){t-/ + t + 277 ) 7 ^ 

(l + t77)3(l-t2) ^ ’ 

6(7 - l)(t7 + t +277)7'^ 

(l + t 77 ) 7 i + V 7 ) 7 l-t")' 

Hence near (ao,ro), the power series of H{a,r) has the form 

H{a,r) = Hra{ag,rg){r - rg){a - ag) + ^Haa{ag,rg){a - Oo)^ + ^Hrrr{ag,rg){r - rg)^ H- 

Motivated by this expansion, we set 


Tdra{.CLg, rg) 

Haa{ag,rg) 


(5.34) 

(5.35) 


(5.36) 


r -rg 

V ■= -> 

(1 — (2q 

and substitute r = rg + rj{a — ag) in H, and define 


(a - ao)2 


(5.37) 


(5.38) 
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Setting 


2+2 


Haa{aQ,rQ) 3 ( 7 -l)(t 7 + t + 2 ^) 2 t' 


2Hra{ao,ro) 4{1 + t^){t + ’ 

a direct calculation shows that F{a,r]) is a smooth function near ( 00 , 70 ) and 

F{ao,Vo) = 0: Fvi^o,Vo) = Hra{ao,ro) ^ 0. 


(5.39) 


(5.40) 


Therefore by the implicit function theorem, there is the smooth function 7 = 7 (a), uq < a < ao(l + S) 
for some ^ > 0 such that F(a, r]{a)) = 0. Then H{a, tq + 7 (a)(a —oq)) = 0, and by the uniqueness of the 
solution, r = ro + ri(a){a — uq) is smooth in a for oq < a < ao(l + (5). But for a > oq, F[a(a,r(a)) > 0, 
and the smoothness of r = r(a) is elementary. Hence r(a) is a smooth function for a > oq. □ 


Let 


xq := x(ao) = 


V7(l-t^) 


2/0 := 2 /(ao) = 


l-F 


which are obtained by setting a = ao, r = r{ao) = tq in (3.25|), (3.26). Then 


^2 _ 2 _ f-2 

■^0 '0 ^ _ _i 

2rot-i “ ’ 


(5.41) 


(5.42) 


and hence from ( 3.351) , the point ^ is on the negative real line. Thus when a = Og, the two endpoints ^ 
and ^ collapse to the point —rg on the real line. This is an extreme case of the deformation illustrated 
in Figure [l^ . 

Proof of Lemma |5.4 When a = ag, we have r = tq, x = xq, y = yo, and we can direct check from 


( 5.26 ) that A(ag) = 0. We have 

A'(a) = - ylogt + (1 + 7 ) log 2 + logr - 7 log(r + + x) - log(r FtFy) 


, r' 1 r' — X 

+ (l + a- 

r 


1 r'+ x' 

2 r + — X 2 r + + x 2 r + t — y 


r' +y' 


^ r-y' W, „ , 

--(l + 2a)-. 

2 ^ r + t + y 


At a = og, from from ( ^.29 ), ( 5.41 ), we have 

- 7 logt + ( 1 + 7 ) log 2 + logrg - 7 log(rg + a;g) - log(rg + t + yo) = 0 

and hence again from ( 5.29 ), ( ^.41 ), after some algebra, 

(7-l)(l-t^) 


A'(ag) — —- 


{t + tj + 2^)il + t^) 

Now from the relation ( tl.27| ) between r, x, y, we have 

1 


[(1 + ty^yiao) + ty^x'iao) - //'(og)]. 


1 -F 


{yy' - t^xx'). 


(5.43) 


(5.44) 


(5.45) 


(5.46) 
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This implies by ( 5.29| ), ( 5.41 ), 


(1 + ty^y{ao) = y\ao) - t^fx^ao), 


and hence 


A'(ao) =0. 


(5.47) 


(5.48) 


Now we compute A"(ao). We have 


A"{a) —2— + (1 + a) ( —— (—) I - 


r\ 2 \ If r"-x'' 




-27 


— -(1 + 27 a) 


r + + X 2 

'! n.f 


2\r + t~^ — X ^r + t~^—x' 
r'' + x'' , r' + x' .2 

— llrr-/ 


r-\-t -*-+ X ^r + t -*-+x^ 




r" + y" 
r + t + y + t + y 


y”y' y) -2 ''' --{i + 2a) 

2\r + t — y ^r + t — y' J r + t + y 2' 


First consider the terms with double derivatives. At a = ag, 

,r" 1 r"-x" 1, „ , r'' + x'' I r” - y” 1,^ r. ^ r" + y" 

(1 + a)---—--(l + 27 a) -(1 + 2a)- 


r 2 r + t ^ — X 2 

h-m-y 


4{t + tj + 2yy){l + tyy) 


r + t ^+a; 2r + t — y 2 
[(1 + tyiViag) + tyjx"{ag) - y”{ao)]. 


r + t + y 


From (5.46), 


rr" - Y^iyy" - = -fr'f + 


(5.49) 


(oo) 


(5.50) 


(5.51) 


and hence the right-hand side of (5.50D is equal to 


- >+’Mf)]. (« 2 , 


Thus A"(ao) is given by (B at a = Co where the terms with double derivatives are replaced by ( 15^^) , 
which involves only the first derivatives of r,x^y at oq. From ([5.30D and (3.25), (p.26), we have 


x'{ao) = 

y'(oo) = 

and we obtain, after some calculation. 


3(7- + t + 2yyff^ 

4(l + tV7)(i + V7)3 

(5.53) 

it + 4/7 -b 3<7)(t + t-i + 2yj)H 

(5.54) 

4(l + tV7)(t + V7)3 

(4tV7 + 7 + 3)(i + ^7 + 

(5.55) 

4(l + tV7)(i + V7)3 

f^it + t-i + 2yyfy^ 

2{i + t^f{t + y^f 

(5.56) 
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By Taylor’s formula, for a > uq, 

A(a) = iA"(ao)(a - ao)^ + - ao)^ (5.57) 

for some ag < a < a. For 5 > 0 and 1 < 72 < 00 , set 

C := sup{|A"'(a, 7 )| : 1 < 7 < 72 , 00 ( 7 ) < o < 00 ( 7 ))! + i5)}, (5.58) 


where we have made the dependence on 7 of A'" explicit. It follows from the smooth dependence of 
A(a, 7) on 7, as well as on a, that given 72, we can choose 6 = S 2 such that C < 00 . Therefore 

A(a) = iA"(ao)(a - ao)^(l + 0{\a - oo|), (5.59) 

where 0(|a —oqI) is uniform for 1 < 7 < 72 and ag < a < ao(l + (52). Indeed 0(|a —oqI) < Clo —ao|. □ 


In order to prove Proposition 2.2 


we use (^.18[) . As 

x{a)‘^ — r(a)^ — 


cos 9c = 


2 r{a)t 


-1 


(5.60) 


and cos 0c ^ — 1 as o —> oq, we see that sin ^ is a smooth function of a in [oq, 00). Thus 

sin ^ = 1 + 0(|a — ttol) for a near ag, a > ag. Inserting this information into ( 5.18| ) and using ( |5.27| ), 
we obtain (^.2C|). 


6 Proof of Theorem |1.1| 

Take 


L > 


2Lo 
4/3 T ’ 

% bg 


0 < s < 


agbg{l + <5o) ’ 


( 6 . 1 ) 


where Lg^Sg are given in Proposition 2.2. Let 


Set 


—N-xbgN^/^. 

(6.2) 

ag 

00 

(6.3) 


Then for L < a: < the condition ( 2.19| ) for a = ^ in Proposition |2.2| is satisfied for any k in the 

range n < k < b, Following |Q, we consider 

,P(G([7A], A) < n) 


log 


P(G([7 A],A)< 6 ) 


which equals 


^ log(-y2i(0;fc)) 


(6.4) 


(6.5) 
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by (^). Inserting ( 2.20| ) into (^), we obtain (cf. |I^ ) 


ImjNimlt = +o(i„g x) 

= + 0{logx) 

as C2aQ6g = j. But by the result of Johansson Q, as N ^ oo, 

P(G([7iV], fV) <b) = Fi-Loa-^/%^) + o(l), 
where F{x) is the Tracy-Widom distribution. This proves Theorem 0 


( 6 . 6 ) 


(6.7) 
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